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Abstract. One possible scenario for the origin of the solar tachocline, known as the "fast tachocline", assumes 
that the turbulent diffusivity exceeds r\ J; 10 9 cm 2 /s. In this case the dynamics will be governed by the dynamo- 
generated oscillatory magnetic field on relatively short timescales. Here, for the first time, we present detailed 
numerical models for the fast solar tachocline with all components of the magnetic field calculated explicitly, 
assuming axial symmetry and a constant turbulent diffusivity r\ and viscosity v. We find that a sufficiently strong 
oscillatory poloidal field with dipolar latitude dependence at the tachocline-convective zone boundary is able to 
confine the tachocline. Exploring the three-dimensional parameter space defined by the viscosity in the range 
log v = 9-11, the magnetic Prandtl number in the range Pr m = 0.1-10, and the meridional flow amplitude (—3 to 
+3cm/s), we also find that the confining field strength B con f, necessary to reproduce the observed thickness of the 
tachocline, increases with viscosity u, with magnetic Prandtl number v/rj, and with equatorward meridional flow 
speed. Nevertheless, the resulting -B con f values remain quite reasonable, in the range 10 3 -10 4 G, for all parameter 
combinations considered here. The thickness of the tachocline shows a marked dependence on both time and 
latitude. The latitude dependence is similar to that inferred by helioseismology, while the time dependence is 
within the observational errors. 
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1. Introduction 

The tachocline is the thin transitional layer below the solar 
1 convective zone between the surface-like differential rota- 
tion pervading the convective zone and the rigid rotation 
of the radiative interior. The existence and properties of 
this layer have bee n known from heliose ismic studies, re- 
cently reviewed bv lCorbard et all (l200ll) . The tachocline 
is known to be extremely t hin jcorbardet allll998lll999t 
ICharbonneau et all ll999). Precise values for its thickness 
depend on the particular fitting profiles used in helioseis- 
mic forward modelling for the residual rotation rate ft — Qq 
(where flo is the rotation rate of the solar interior). The 
scale height of f2 — Slo is generally found to be crudely 
0.01-0.02 Rq, so that the residual rotation is reduced by 
one order of magnitude within a layer of barely ~ 30 Mm. 
Near the equator this layer lies just below the adiabatically 
stratified convective zone. Inversions seem to suggest that 
at higher latitudes it is situated at a slightly but signifi- 
cantly higher level, and is apparently also thicker, so that 
it par tly overlaps with the convective zone l|Basu fe AntiaL 
1200 lh . 



The extreme thinness of the tachocline implies a 
strongly anisotropic (predominantly horizontal) transport 
of angular momentum. Several different mechanisms have 
been proposed for this transport, but it is now widely be- 
lieved that the magnetic field is instrumental in its origin. 
Magnetic fields can lead to the required anisotropic trans- 
port either directly, via the action of Maxwell stresses, 
or indirectly, by rendering the rotation profile unsta- 
ble, a nd thus giving rise to an i sotropic turbu lent trans- 
port l|Gilman fc Dikpati 1200(1 iGilmanl |2000|) . (In con- 
trast, purely hydrodynamic mechanisms are apparently 
unable to provide the necessary l a teral transport — cf. 
ICharbonneau fc MacGregorl l|l997|) : iGaraudl (|2001a|U 

Depending on the value of the magnetic diffusivity, this 
magnetic field may cither be a weak permanent, primor- 
dial field pervading the solar interior, or the strong oscil- 
latory field generated by the solar dynamo. A magnetic 
field oscillating with a circular frequency w cyc = 2n/P, 
P = 22 years is known to penetrate a conductive medium 
only down to a skin depth of 
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where rj is the magnetic diffusivity. Comparing this with 
the tachocline scale height quoted above, it follows that 
for 77 <J 10 8 cm 2 /s the dynamo field cannot penetrate the 
tachocline, and we can expect the tachocline to be per- 
vaded by the internal primordial field. On the other hand, 
for j] <; 10 9 cm 2 /s the tachocline dynamics should be gov- 
erned by the dynamo field. As the associated diffusive 
and Lorentz timescales are also very different, these two 
cases basically correspond to the case of "slow" and "fast" 
tachocline, discussed in th e literature (see esp. Table I in 
IGilmanL EoOfl IBrunl Eool . 

The case of a slow tachocline llGoueh fe Mclntvrel 

1998|) has been studied extensively in recent years 
by a number of authors. In models w it h no 
meridional flow llMacGregor fc Charbonneail [l999: 
iRudiger fc Kitchatinovl [l997) it was found that for 
purely molecular diffusivities even a rather weak pre- 
scribed internal poloidal field of G is sufficient to 
confine the tachocline to its observed thickness if the 
internal field is fully contained within the radiative zone. 
The steady toroidal field resulting from the winding up 
of poloidal field lines may reach kilogauss values. The 
containment of the internal field is, however, not so 
easy to explain, given that the surrounding convective 
zone, with its high diffusivity, is far from being a perfect 
conductor. On the other hand, for an internal field not 
contained within the radiative zone Ferraro's law of 
isorotation predicts that the differential rotation will 
penetrate deep into the radiative zone. This wa s indeed 
confirm ed in the numerical models. Recently, lOaraudl 
I 2001b) performed calculations taking into account the 
meridional flow and the self-consistent evolution of the 
poloidal field (with otherwise simplified physics). Her 
models show that there exists an intermediate range of 
field strengths where the interaction with the meridional 
flow indeed leads to the poloidal field being mostly 
confined to the radiative zone, with a resulting rotation 
profile that is comparable to the observations. At any 
rate, these calculations stress the importance of including 
meridional flows and a self-consistent treatment of the 
poloidal field. 

The alternative case of a fast tachocline, in con- 
trast, has received much less attention. Th is is so de- 
spite that several authors, fe.g. IGilmanL 12000^) noted that 
the tachocline is more likely to be dynamically coupled 
to the convective zone than to the radiative zone, and 
that a part of the tachocline is likely to overlap with the 
strongly turbulent convective zone, especially at higher 
latitudes. Faced with this situat ion, in a recent paper 
l|Forgacs-Daika fc Petrovavl (|200l[) . hereafter FDP01) we 
considered the case of a turbulent tachocline with a tur- 
bulent diffusivity of about r\ = 10 10 cm 2 /s, pervaded by 
an oscillatory dipole field. This model can be regarded 
as the anal ogue of the models of IRudiger fc Kitchatinovl 
l)l997|) and iMacGregor fc Charbonneaul l)l999l) for a fast 
tachocline. It was found that for a sufficiently high value of 
the field strength, Maxwell stresses in the oscillatory field 
transported angular momentum efficiently enough to con- 



fine the tachocline to its observed thickness. The poloidal 
field strength necessary for this confinement could also be 
estimated from a simple analytic relation; with the diffu- 
sivity quoted above its value was 2400 G. 

The analysis of FDP01 was of limited scope, as sev- 
eral simplifying assumptions were made. Only the case of 
magnetic Prandtl number Pr m = v jr\ = 1 (v is the viscos- 
ity) was considered, meridional circulation was neglected, 
the poloidal field was prescribed throughout the computa- 
tional domain, the applied field was an oscillating dipole 
instead of a dynamo wave, the diffusivities were arbitrarily 
prescribed as constant or known functions of radius etc. 
The aim of the present series of papers is to systematically 
generalize the model of the fast tachocline by incorporat- 
ing more and more of the effects neglected in FDP01. As 
a first step, in the present paper we will incorporate the 
meridional flow and include a self-consistent calculation 
of the poloidal field, t he impor t ance of which was under- 
lined by the results of iGaraudl l(2001bh in the context of 
the slow tachocline. We also explore the parameter space a 
bit more widely by varying the magnetic Prandtl number 
and the diffusivity. 

After formulating the mathematical problem in Sect. 2, 
results of the model calculations will be presented in 
Sect. 3. Finally, the Conclusion summarizes the main re- 
sults. 

2. The model 

2.1. Equations 

The time evolution of the velocity field v and magnetic 
field B are governed by the Navier-Stokes and induction 
equations. We write these in a frame rotating with the 
fixed internal rotation rate Qq: 

d t v +(v • V)v = V Q|« x if -V^j - 2(fl x v) 

--V (p+ ^) + J-(B • V)B + -V -f, (1) 

and 

<9 t B = Vx(vxB)-Vx (77V x B), (2) 

where r is the radius vector, f is the viscous stress tensor, 
77 is the (turbulent) magnetic diffusivity, V is the gravi- 
tational potential and p is pressure. These equations are 
supplemented by the constraints of mass and magnetic 
flux conservation 

V • (pv) = 0, (3) 

V ■ B = 0. (4) 

We assume spherical geometry and axisymmetry in all 
our calculations. The velocity and magnetic fields can be 
written as 

v = uj(r,9,t)r sin Oe^ + v m (5) 
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B = B(r, 8, t)e^ + V x [A(r, 8, t)e 4 



(6) 



where B(r,9,t) and A(r, 9,t) represent the toroidal com- 
ponents of the magnetic field and of the vector potential, 
respectively; u)(r, 9,t) = O — il is the residual angular ve- 
locity and v m = v r e r + vgeg is the meridional circulation. 

As a self-consistent calculation of the meridional circu- 
lation throughout the dynamically coupled regions of the 
convective zone and the tachocline is beyond the scope 
of this work, we will regard the meridional circulation as 
given. As a result of the anelastic approximation, the cir- 
culation in the spherical shell may be represented by a 
stream function so that the components of the merid- 
ional circulation can be written as 

Vr = -^—.^ (7) 

pr* sin 9 

ve = —^d,^. (8) 
pr sin 9 

In order to present more transparent equations we write 
the stream function in the following form: 

* = i/j(r) sin 2 6 cos 6. (9) 

the form of the given function ip{r) specifying the flow. 

The evolution of the poloidal magnetic field is, in con- 
trast, self-consistcntly calculated by integrating Eq.0 as 
usual: 

dtA +- 1 r- ^—--a^--./^ 1 



-(v m -V)(rsin0A)=77 V 2 - 



r sin 9 



A.(W) 



The components of the viscous stress tensor appearing 
in the azimuthal component of the Navier-Stokes equation 
are 



T4>e = pv- 



T rt j, = pur 



r Vsi 



sin 8 J ' 



(11) 
(12) 



where v is the (turbulent) viscosity. With this, the az- 
imuthal components of Eqs.EEl including the effects of 
diffusion, Coriolis force, meridional circulation, toroidal 
field production by shear, and Lorentz force, read 



d t uj 



o , a v , dr P\ o , Zvcote 
a r v + 4 — V v ) o r uj H s — dguj 



(13) 



+ vd*u + ^d 2 e u + L + Mi + C 



L 



B 



<iirpr 2 sin 9 

d r B 
Airpr 2 sin9 



d e A 



cot 9 d r A 



Airpr 2 sint 



{d e A-cot9A) 



- + d r A 



1 ( ~ 9 „ n , 1-3 cos 2 6 „ , , 

Mi = — 2 cos 2 8 d r ip H 2ip ui 

r l p \ r 



1-3 cos 2 i 
r 2 p 



ip d r u> 



sin a cos ( 
r 2 p 



d r ip dgu 



r m ° n q 2* , , 2» cos 2 g Q 
C = — =— (1-3 cos 6*) V H o d r ip 



d t B = I — 



M 2 = 



r 2 sin 2 i 



B 



— +d r r 1 \ d r B 

r 



(14) 



V 



n 



cot9 d 9 B + rjd 2 B + 8 2 B + M 2 



(3cos 2 0-l) I - + — 
1 r p 



B 



r 2 p 

_ <W cos 2 9B + ^L(l-3cos 2 9)d r B 
r 2 p r 2 p 

while equation (|lt)|l can be written as 



d t A = 



-V 



r 2 sin 2 i 



A 



2q 



d r A 



r\ cot 8 



d 9 A 



(15) 



•Ma 



1 / 1 - 3 cos 2 8 



2 



+ cos 2 8 d r ip A 



/i o _2 - ■ - sin cos ( 



-^(1 - 3 cos 2 6>)<9 r A 



r 2 p 



d r ip dgA 



In these formulae L denotes the terms associated with the 
Lorentz force, C denotes the terms associated with the 
Coriolis force and Mi, M 2 , M3 denotes the terms associ- 
ated with the advection by meridional circulation. 

2.2. Boundary and initial conditions 

The computational domain for the present calculations 
consists of just the upper part of the radiative interior, 
between radii n„ and rb cz , where rb cz is the radius of 
the bottom of the convection zone. For the integration of 
Eas. ll3tfl5l we use the following boundary conditions. 

At the pole and the equator axial and equatorial sym- 
metry is required for w, and dipole symmetry for the mag- 
netic field: 



d e oJ = B= A = 
d H uj = B = d H A = 



at 9 = 0, 

at 8 = tt/2. (16) 



The bottom of our box is supposed to be a rigidly 
rotating perfect conductor — a good approximation for 
the solar interior. The rigid rotation speed is fixed at the 
observed rate Vt /2ir — 437 nHz. Thus, the lower boundary 
conditions are: 



lu = B = A = 



at r 



(17) 



Note that the fixed internal rotation rate Oo is in general 
not equal to the rate at which the interior of our mod- 
els would rotate if the lower rotation rate were not pre- 
scribed. (In other words, the fixed rate does not in general 
ensure zero net torque at the lower boundary.) This leads 
to a spurious radial differential rotation in the radiative 
zone below the tachocline. In models where this effect is 
strong, it is a sign that the model does not represent well 
the observations as it would spontaneously lead to an in- 
ternal rotation rate that differs significantly from what is 
observed. As, however, the magnetic field geometry con- 
sidered here is still not realistic, firm conclusions on this 
would be premature anyway. 
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The lower boundary needs also to be placed sufficiently 
deep in the solar interior to make sure that the confine- 
ment of the tachocline is not an artefact caused by the 
rigid rotation imposed at the bottom. For this reason, in 
most calculations presented here we set r in = O.4R , even 
though only the upper part of our computational domain 
is shown in many of the figures. 

At the upper boundary of our domain, i.e. at the bot- 
tom of the convective zone, we use the following boundary 
conditions: 



— ^bcz — 

A = Aq sin 9 cos(wcyct) 
B = 



at r — Tbczj 



(18) 



The physical assumptions behind these conditions are the 
following. For u> we suppose that the rotation rate can be 
described with the same expression as in the upper part of 
the convection zone. In accordance with the observations 
of the GONG network, the expression used for f^bcz is 



Ob 



2tt 



= 456 - 72 cos 2 6>-42 cos 4 6 nHz. 



(19) 



Our upper boundary condition on the poloidal field is the 
simplest possible representation of the poloidal compo- 
nent of the solar dynamo field, a dipole oscillating with a 
period of 22 years. The factor Aq fixes the poloidal field 
amplitude. Finally, the upper boundary condition chosen 
for the toroidal field assumes that the toroidal field is neg- 
ligible in the convective zone compared to its value in the 
tachocline. Such a hypothetical situation is in line with 
current thinking in dynamo theory, and it may be a natu- 
ral consequence of buoyancy-driven instabilities effectively 
removing any toroidal flux from the convective zone. 
The initial conditions chosen for all calculations are 



(20) 



Some test runs with other initial conditions (with signifi- 
cant differential rotation rotation throughout the volume) 
were also performed. After stationary periodic behaviour 
sets in, the results were found to be independent of the 
initial conditions chosen. 



2.3. Numerical method 

We used a time relaxation method with a finite difference 
scheme first order accurate in time to solve the equations. 
A uniformly spaced grid is set up with 128 grid points in 
the r direction and 32 grid points in the 9 direction. 

Our calculatio ns are based on a more recent version of 
the solar model of iGuenther et al" 

Starting from the initial conditions 1)20(1. the solution is 
allowed to evolve in time until it relaxes to a very nearly 





9,t 


= 0) 


— Obcz — ^0 


at T — Hjcz 


uj(r,( 


9,t 


= 0) 
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at r < r hcz 


B{r, 


9,t 


= 0) 


= 




A(r, 


9,t 


-o) 


= A sin 2 9 


at r = rbcz 


A(r, 


e.t 


= 0) 


= 


at r < r bcz 



periodic behaviour. The timescale needed for this is the 
diffusive timescale 



idif 



(?"bcz ^*in) 

Min(i/, 7]) 



(21) 



where Min(i^, rj) is the lowest value of viscosity or mag- 
netic diffusivity in the domain. Physically we would expect 
this value to be close to the molecular viscosity. Using 
this value in the computations would, however, lead to 
a prohibitively high number of timesteps to relaxation. 
(Clearly, the runtime of our computations must be chosen 
to well exceed idif to reach relaxation.) Consequently, in 
the present paper we only consider cases where the diffu- 
sivities are constant throughout the domain. Note, how- 
ever, that in FDP01 a case with the diffusivity significantly 
reduced towards the interior was also computed, and no 
qualitative changes in the character of the solution were 
found. 

The question may arise, to what extent does the use 
of such a high viscosity in the deep radiative zone distort 
the results? The answer is clearly that the main differ- 
ence between models with high and low diffusivities is the 
diffusive timescale over which the model evolves to equi- 
librium. Once equilibrium is reached, the azimuthal com- 
ponent of the equation of motion reduces to a Laplace 
equation (neglecting meridional circulation and magnetic 
fields), so the solution is universal and independent of the 
diffusivity. The equilibrium solution is thus not expected 
to be strongly influenced. 

On the other hand, for realistically low (such as molec- 
ular) diffusivities equilibration takes much longer than the 
age of the Sun. Thus, an internal differential rotation, in- 
duced by solar wind torques or as a pre-MS relict, may 
not be suppressed. It is therefore important to make it 
clear (as we already did in FDP01) that the present work 
does not propose that dynamo-generated fields can ex- 
plain the overall lack of differential rotation in the deep 
radiative zone, just that they dominate the dynamics of 
the tachocline itself. The lack of differential rotation in 
the deep radiative zone still needs to be explained by an 
internal magnetic field independent of the dynamo. 

3. Results and discussion 

3.1. Reference case 

First we present results for a case with Pr m = 1, rj = 
v = 10 10 cm 2 /s, and no meridional circulation (-0 = 0). 
This "reference case" corresponds to that computed in 
FDP01, except that the poloidal field is now not pre- 
scribed throughout the volume, but explicitly calculated. 
The results are shown in Fig.^ 

In the left-hand panels we plotted the contours of the 
time-average of the angular rotation rate, which is defined 
as 



1 



Pi 



eye Jt 



t+Pcyc 



w(r, 0,t)dt. 



(22) 
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Fig. 1. Spreading of the differential rotation into the radiative interior for r\ — v = 10 cm /s. Left-hand panel: 
contours of the time-average of the angular rotation rate u(r,8,t) under one dynamo period. Equidistant contour 
levels are shown, separated by intervals of 100nHz/f2o, starting from towards both non-negative (solid) and negative 
(dashed) values. Right-hand panel: normalized differential rotation amplitude Aw at different latitudes as a function 
of radius. The peak amplitude of the poloidal magnetic field is B p ~ 2570 G. 
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Fig. 2. Same as in Fig.^for rj = v = 10 11 cm 2 /s. The peak amplitude of the poloidal magnetic field is B p ~ 6660 G. 
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Fig. 3. Same as in Fig.^for r\ = v = 10 9 cm 2 /s. The peak amplitude of the poloidal magnetic field is B p ~ 1386 G. 
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Fig. 4. Snapshots of the solution at four cycle phases for r) = v = 10 10 cm 2 /s. Left-hand panels: contours of poloidal 
magnetic field strength. Equidistant contour levels are shown, separated by intervals of 500 G. Middle panels: contours 
of toroidal magnetic field strength. Equidistant contour levels are shown, separated by intervals of 1000 G. Right-hand 
panels: contours of angular rotation rate, as in Fig.Q] 
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In the right-hand panels the differential rotation am- 
plitude Aid is defined as 



Aw = w(r,6,t) 



Aw = w(r, 9) - [ 
Jo 



r/2 



u(r,9,t) sin0 d6 



r/2 



w(r,f 



) sin 9 d9. 



(23) 



(24) 



The thickness of the tachocline will be defined throughout 
this paper as the scale height of Aw (i.e. the depth from 
the top of our domain where Aw is reduced by 1/e). By 
"mean thickness" , in turn, we will refer to the scale height 
of Aw. 

The results shown in Fig.^ are in accordance with 
those in FDP01. In particular, we confirm that an oscilla- 
tory poloidal field with dipolar latitude dependence at the 
tachocline-convective zone boundary is able to confine the 
thickness of tachocline. The effectivity of this confinement 
increases (i.e. the thickness of the tachocline decreases) 
with the field strength. 

This conclusion is basically valid for all the cases stud- 
ied in this paper, though the field strength necessary to 
just reproduce the observed equatorial mean thickness of 
the tachocline, to be called the confining field strength 
-B con f, for brevity, varies from case to case, as do the struc- 
tural details of the resulting tachocline. Though in or- 
der to find the right value of -B C onf we experimented with 
many different field strengths for each parameter combi- 
nation, in this paper we only present results computed 
with B = B con f (cf. FDP01 for calculations with other 
field strengths). The value of the confining field strength 
is found to be ~ 2600 G for the reference case, very close 
to the value found in FDP01. 

In what follows, we will systematically study the effects 
of varying the parameters of the model on the resulting 
tachocline structure and on the value of the confining field 
strength. First, we fix the value of the magnetic Prandtl 
number as Pr m = vjr\ = 1.0, while we change the ampli- 
tude of the magnetic diffusivity. Second, we examine the 
influence of varying the magnetic Prandtl number on the 
radial spreading of the differential rotation into the radia- 
tive interior. In these two cases we ignore the meridional 
flow. Finally, we will focus on the effect of introducing a 
meridional flow. 

Computer animations illustrating the time develop- 
ment of some of our solutions can be downloaded from 
the following web site: 

http : / / astro . elte .hu/kutat/ sol/fast 1/fastle .html 
3.2. Varying the difFusivity 

In this subsection, the value of the magnetic Prandtl num- 
ber was fixed, while we studied the effects of varying the 
diffusivity. The results are shown in Figs. lll3l 

In line with the analytic estimates of FDP01, we find 
that -Bconf increases with the value of the diffusivity. Given 
that the equipartition field strength in the deep convective 
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Fig. 5. The thickness of the tachocline at different lati- 
tudes as a function of time for the case in Fig.^ 
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Fig. 6. Latitudinal variation of the mean tachocline thick- 
ness for the case in Fig.^ A vertical dashed line indicates 
the poleward limit of reliable helioseismic data. 



, :on f is quite reasonable even 
f , a higher diffusivity 



zone is ~ 10 4 G, the value of B. 
with a diffusivity of 10 11 cm 2 /s. 

Beside increasing the value of B con 
has other effects on the solution, too. The fast tachocline 
shows a quite marked time- and latitude-dependence in 
our models, as illustrated in Fig.^for the reference case. 
As a consequence, the thickness of the tachocline de- 
pends both on latitude and cycle phase. In particular, the 
tachocline tends to be significantly wider at higher lati- 
tudes on the mean, while it is also more time-dependent 
near the poles (Fig.EJ. While the relative shortness of seis- 
mic data sets does not permit firm conclusions about the 
cycle dependence of tachocline properties yet, the overall 
form of the mean thickness-latitude relation (Fig.jBJ is in 
good qualitative agreem ent with the ava ilable helioseismic 
information (cf. Fig. 9 in lBasu fc Antiall200ll) . The ampli- 



8 



E. Forgacs-Dajka and K. Petrovay: Dynamics of the fast solar tachocline 



Pr m =10.0 
T] = 1 0'cm/s 




1 .0 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 
8 [rodl 




Fig. 7. Same as in Fig.Q]for Pr m = 10.0, where r) — 10 9 cm 2 /s and v = 10 10 cm 2 /s. The peak amplitude of the poloidal 
magnetic field is B p ~ 8531 G. 
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Fig. 8. Same as in Fig.^for Pr m = 0.1, where r\ = 10 10 cm 2 /s and v = 10 9 cm 2 /s. The peak amplitude of the poloidal 
magnetic field is B p ~ 1028 G. 



tude of these variations within the latitude interval where 
the helioseismic results are reliable (9 ^ 20°) also seems 
to be in agreement with the observations. 

For a fixed value of the field strength, a higher value of 
diffusivity tends to reduce the amplitude of the spatiotem- 
poral variations. On the other hand, for higher diffusivities 
Bconf is also higher, so for this field strength a stronger dy- 
namical variation is imposed on the flow. Nevertheless, at 
least for the reference case, the time dependence is well 
within the observati onal errors as shown e.g. in Fig. 4 of 
iBasu fc Antial l)200l|) . 

It must be noted of course that, as our imposed upper 
boundary condition on the poloidal field, a simple oscillat- 
ing dipole, may not be very realistic, the detailed spatio- 
temporal properties of the solution should be regarded 
with caution. Among such features we may mention the 
poleward drift of the resulting rotational pattern, appar- 
ent e.g. in the phase lag between different latitudes in 
Fig.El (Cf. the same trend in Fig. 10 of FDP01.) 



3.3. Varying the magnetic Prandtl number 

Figs.[3and[5]illustrate the effect of magnetic Prandtl num- 
bers other than 1 on the solution. As expected on the basis 
of the analytical estimates of FDP01, tachocline confine- 
ment is less efficient for high magnetic Prandtl numbers, 
so that higher values of -B C onf are necessary if Pr m > 1. 
On the other hand, a higher magnetic Prandtl number 
(i.e. a relatively higher viscosity) is in general beneficial 
for the overall smoothness of the solution, leading to less 
strong dependence of uj on time and latitude. In the low 
Pr m case, in turn, strong spatiotemporal variations in to 
become prevalent, leading to a complex non-monotonic be- 
haviour of the rotational velocity as a function of radius 
and latitude. 



3.4. Influence of the meridional flow 

In the case of a fast tachocline dynamically coupled to 
the convective zone, as assumed here, we expect that any 
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Fig. 9. Streamlines of the meridional circulation pre- 
scribed in Eq. Q25[l. with the following parameter values: 
V> = -10 20 ,/c= l,ro = (i? -r mc )/3O andT = 3.4xl0 10 . 
Dotted lines represent counterclockwise circulation. The 
amplitude of the flow is ~ 3 cm/s in our computational 
domain (between the dashed and dotted circles). 




Fig. 10. Same as in Fig.Hwith ip = 10 20 , k = 2. Solid 
and dotted lines represent clockwise and counterclockwise 
circulation, respectively. 



meridional circulation present in the tachocline is unsepa- 
rable from the complex problem of meridional circulation 
in the convective zone. It is far beyond the scope of this 
paper to attempt to model this circulation in a consistent 
way. Instead, we will limit ourselves to regard the effects 
of a meridional circulation with a simple and plausible, 
prescribed spatial structure, and a realistic amplitude. 

It is known from direct and seismic studies that the 
meridional circulation near the phot osphere is pol e ward, 
with a peak am plitude of ~ 20m/s ijKomm et all Fl993: 
iLatushkoL Il994h . If we assume that the meridional flow 
in the convective zone+tachocline system follows a simple 



one-cell pattern, the circulation in the tachocline will be 
equatorward. The alternative possibility of a two-cell pat- 
tern with poleward flow i n the tachocline, how ever, should 
also be left open for now l)Kiiker fc StbJ . l200ll) . An expres- 
sion for the function ip(r), defined in Eq.El that encom- 
passes both possibilities is 



ip(r) = -00 sin 



kw(r 



exp 



- Tgf_ 

r 2 



(25) 



where ipo sets the amplitude of the meridional circulation, 
ro , r and k are geometric parameters and r mc is the radius 
value down to which the meridional flow penetrates from 
the base of the convective zone. Keeping in mind the ob- 
served depth of tachocline, here we set r n 



4.6xl0 10 cm 



The strong subadiabatic stratification of the solar inte- 
rior sets an upper limit to the amplitude of the meridional 
flow. The timescale of any meridional circulation cannot 
be shorter than the relevant thermal diffusive timescale 
(to allow moving fluid elements to get rid of their buoy- 
ancy). In the present case the relevant heat conductiv- 
ity is the turbulent one, so for turbulent Prandtl num- 
bers not very different from unity, the shortest possible 
timescale for meridional circulation is just comparable 
to the viscous timescale. An elementary estimate yields 
v max ~ r bc zAdif ~ 10 cm/s. 

The flow parameters used in our calculations were 
chosen to produce simple smooth one- and two-celled 
flow patterns obeying the amplitude constraint discussed 
above, while reproducing the observed flow speed near the 
surface. The flow patterns are shown in Figs.l9landllt)l In 
these cases, the speed of the horizontal motion is ~ 3 cm/s 
in the upper part of the radiative interior. 

The results with the meridional circulation included 
can be seen in Figs.^2 an d 03 The poloidal field ampli- 
tude used here was the same as in the reference case. It 
is apparent that an equatorward flow in the tachocline re- 
gion reduces somewhat the efficiency of the confinement 
(the tachocline becomes slightly thicker), though the ef- 
fect is not dramatic. A poleward flow, on the other hand, 
significantly improves the efficiency of the confinement. 

In fact, meridional circulation alone, without any im- 
posed magnetic field, can also have a significant effect on 
the thickness of the tachocline, as it can be seen by com- 
paring Figs.1131 and 1141 

The relative importance of meridional circulation 
scales approximately with the ratio of the diffu- 
sive/advective timescales vjv\ thus, for a higher value of v, 
a proportionally higher flow speed is needed to reproduce 
the same tachocline thickness. 

4. Conclusion 

We have begun a systematic study of the dynamics of 
the fast solar tachocline arising under the assumption 
that the turbulent diffusivity in the tachocline exceeds 
i] il 10 9 cm 2 /s. We extended the earlier preliminary anal- 
ysis of FDP01 by consistently computing the evolution of 
the poloidal field, and by exploring the three-dimensional 



10 



E. Forgacs-Dajka and K. Petrovay: Dynamics of the fast solar tachocline 



0.65| 

0.60 
°? 0.55 
0.50 
0.45 



Pr m =1.0 
7) = 1 0'°cm/s 
counterclockwise 1.0 
circulation 



y 



0.07 
0.01 
-0.05 

c 

\ 

-0.1 1 

0.17 
0.23 
0.29 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 
8 [radl 




Fig. 11. Same as in Fig.^ but including the meridional circulation shown in Fig.|5J 
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Fig. 12. Same as in Fig.^ but including the meridional circulation shown in Fig.llQI 



parameter space defined by the viscosity in the range 
logjy = 9-11, the magnetic Prandtl number in the range 
Pr m = 0.1-10, and the meridional flow amplitude (—3 to 
+3cm/s). The feedback of the Lorentz force to the merid- 
ional circulation was neglected. 

Confirming the results in FDP01, we found that, for 
basically all parameter combinations considered, an oscil- 
latory poloidal field with dipolar latitude dependence at 
the tachocline-convective zone boundary is able to confine 
the thickness of tachocline. (Note, however, that for low 
magnetic Prandtl numbers the flow becomes rather com- 
plex, and a simple, smooth, monotonic tachocline layer 
does not form.) The effectivity of this confinement in- 
creases (i.e. the thickness of the tachocline decreases) with 
the field strength. The confining field strength B con { neces- 
sary to just reproduce the observed mean equatorial thick- 
ness of the tachocline is found to increase with increas- 
ing viscosity, magnetic Prandtl number v/rj, and merid- 
ional flow speed (if considered positive equatorwards). 
Nevertheless, the resulting B con { values remain quite rea- 
sonable, in the range 10 3 -10 4 G, for all parameter combi- 
nations considered here. 



The detailed dynamics of the solutions indicates that 
the thickness of the tachocline, in general, increases with 
latitude, and it also shows a rather strong dependence on 
cycle phase (Figs.|3] and EJ). The latitude dependence is 
similar to that inferred by helioseismology, while the time 
dependence is within the observational errors. 

The models with no meridional circulation also show 
a tendency for the rotational modulation to drift pole- 
wards (cf. Fig.[SJ). It must be noted, however, that owing 
to the rather simplified geometry of the poloidal mag- 
netic field imposed on the upper boundary (oscillating 
dipole), the relevance of our results to the real Sun is 
dubious, as far as the detailed spatiotemporal structure 
of the tachocline is concerned. Another example of such 
features is the polar "pit" of slow rotation in nearly all 
of our models: again, the reality of this pit needs to be 
verified with other poloidal field geometries. A next step 
in our study of the fast tachocline should clearly be an 
extension of this investigation to more general magnetic 
field structures, reminiscent of the butterfly diagram. 

Another important shortfall of these models is obvi- 
ously the extremely simplified treatment of the turbulent 
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Fig. 13. Spreading of the differential rotation into the radiative interior including the meridional circulation shown in 
Fig-EI but taking the imposed field to be null. 




transport, by constant scalar diffusivities. The turbulence 
present in the tachocline is presumably the product of 
complex MHD instabilities. The study of the non linear 
development of the se instabilities has just begun ijCallvl 
l200lUMiesc5l200lj) . and we are clearly a long way from 
being able to give a consistent description of tachocline 
turbulence. Nevertheless, the present models could still 
be significantly improved by considering simple, generic 
cases such as a scalar diffusivity coupled to the value of 
the local shear. Extensions of the present work along these 
lines will be treated in upcoming papers of this series. 
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